FIAN/TD-10/05 
ITEP/TH-45/05 



Matrix Models vs. Matrix Integrals 



A. Mironov a 

Lebedev Physics Institute and ITEP, Moscow, Russia 



• I n a brief review, we discuss interrelations between arbitrary solutions of the loop equations that 

' describe Hermitean one-matrix model and particular (multi-cut) solutions that describe concrete 

matrix integrals. These latter ones enjoy a series of specific properties and, in particular, are 
' described in terms of Seiberg-Witten-Whitham theory. The simplest example of ordinary integral 

r*j , is considered in detail. 
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1 Introduction 

(N ; 

Recent interest in matrix models and especially in their multi-cut solutions was inspired by the studies 
in Af = 1 SUSY gauge theories due to Cachazo, Intrilligator and Vafa and by the proposal of 
Dijkgraaf and Vafa [2] to calculate the low energy superpotentials, using the partition function of 
multi-cut solutions. The solutions themselves are well-known already for a long time (see, e.g., j3]) 
■ with a new vim due to the paper by Bonnet, David and Eynard [1]. 
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The matrix model under consideration is the Hermitean one matrix model. Its partition function 
is given by the integral over N x N Hermitean matrix 

Z N (t) = — ^ / DMexp (tr J> fc M fc ] (1) 

Here DM is the invariant (Haar) measure on Hermitean matrices, which is just the flat measure 
and Volj7(jv) is the volume of the unitary group U(N) 0. Since the integrand in Q is invariant 
w.r.t. matrix conjugation, one can reduce formula (JTJ) to an integral over eigenvalues of the matrix M 
integrating out the angular variables. It can be done using formulas from and the result reads 

Zn{t) = w.I n^o^p (E**** I ( 2 ) 

i \k,i / 

where A(x) is the Vandermonde determinant, A(x) = detjj x{~ 1 = Yli>j( x j ~ x i)- 

The integrals (0) and (J2J) still need to be defined. Indeed, one has to fix the integration contours in 
these multiple integrals. Moreover, these contours can differ from each other. Saying this, one definitely 
forgets that we originally started from Hermitean matrices. However, it is in no way important for the 
properties of Q and ©. In fact, the word "Hermitean" has no other sense but fixing the integration 
contour. What we really need to fix in Q is the flat measure. Say, this is enough to obtain @ from 

©• 

In order to define © and @ one can substitute them with their saddle point approximations 
01 [7] . The other possibility [HJ El El is to observe that they satisfy an infinite set of loop (Virasoro) 
equations (=Schwinger-Dyson equations, =Ward identities) Indeed, this is the consequence 

of only the flatness of measure and invariance of integrals w.r.t. the change of variables. 
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The second possibility means that one calls the matrix model partition function any solution to 
the loop equations. Then, the partition function is not a function but a formal D-module, i.e. the 
entire collection of power series (in t- variables), satisfying a system of consistent linear equations. 
Solution to the equations does not need to be unique, however, an appropriate analytic continuation 
in t-variables transforms one solution to another, and, on a large enough moduli space (of coupling 
constants t), the whole entity can be considered, at least formally, as a single object: this is what 
we call the partition function. Naively different solutions are interpreted as different branches of the 
partition function, associated with different phases of the theory. Further, solutions to the linear 
differential equations can be often represented as integrals (over spectral varieties), but integration 
"contours" remain unspecified: they can be generic chains with complex coefficients (in the case 
of integer coefficients this is often described in terms of monodromies, but in the case of partition 
functions the coefficients are not required to be integer). 

Therefore, in further consideration we distinguish between matrix model, which is a set of solutions 
to the loop equations, and matrix integral which is an integral with a specified integration contour. 
Note that the matrix model partition function is an arbitrary linear combination of a proper set of 
basis matrix integrals. 

The natural choice for this basis is given by the multi-cut (or Dijkgraaf-Vafa) solutions mentioned 
above. They are distinguished by a special property of isomonodromy that allows one to associate 
a Seiberg-Witten-Whitham system with them |131 114j . the corresponding partition function having a 
multi-matrix model integral representation [1J [7] . 



2 lxl matrix case: a toy example 

We start with the simplest example of the lxl matrix, i.e. the matrix integral is just an ordinary 
one-fold integral. Let us start first from the integral with potential of the very general form and see 
what can be done. 



2.1 Loop equations 

Thus, we start from the integral 



Z({t k }) = J^dxexp \ J2t k x k ^J 



(3) 



At the moment we do not specify the integration contour C, but assume it is closed, or, at least, the 
integrand is canceled at the ends of the contour with any monomial x k of arbitrary degree. Then, 
one can immediately obtain an infinite set of equations 1 satisfied by this integral. To obtain it, one 
suffices to consider the integral of the full derivative 



dx-^- 

c ox 



„n+l 



exp tkX k ^J 



n > -1 



which is zero. (@J) gives rise to the set of constraints 

d k+n Z . .dz 

> kt k — r— + (n+l)- 



0, n>-l 



where |^ = Z. Note that 



dZ _ d k Z 
dfk'^tf 



(4) 



(5) 



(6) 



1 These equations are called loop equations, Virasoro constraints, Schwinger-Dyson equations, Ward identities, see 
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and, therefore, the first terms of the sum in © can be rewritten as the linear derivatives w.r.t. the 
couplings t k+n . 

Let us try to solve @ as a power series in t^s, 

Z = c*°> + + (7) 

k k,l 

Then, from the constraint with n=0 we immediately obtain = 0, similarly, cfi = etc. This 
means the solution is trivial: Z = 0. In fact, this should not come as a surprise, since the dimension of 
coupling tk is equal to k, and the dimension of should be negative. However, we have no quantities 
with any negative dimensions at hands. 

In order to get a non-trivial solution, we need to allow, at least some couplings t& to be in the 
denominator. Let us fix a few (say, p) first couplings not to be small, i.e. we shift these couplings 
tk — > + tk, k = l,...,p and consider Z as a power series in t^'s but not in T^'s. Then (JSJ) has 
non-trivial solutions. 

In other words, we fix a polynomial V(x) = Y% ^k xk an d consider the integral in © as taken 
perturbatively w.r.t. to t^'s, 



Z({tk}) = J c dxexp (v(x) + J2hx k j 



(8) 



i.e. one calculates the moments 

' dxx k e v(x) (9) 



c 



The generating function (resolvent) for these moments is 



G(z) = f dx^- = £ f dxx * e V(*) (io) 



c z — x ~ r T1 Jc 



where integral is understood as the principal value integral. 



2.2 "Matrix" model 

Below we consider the simple example of the cubic potential V{x). By the shift of x one can erase 
the quadratic term in the potential and, by rescaling x, make T3 equal to —1/3. Therefore, we are 
left with the only essential variable, which we choose at the moment to be T\. Constraints © taken 
at zero couplings reduce then to the single equation 

g + T lZ = (11) 

which is nothing but the Airy equation. It has two solutions 

Z(T 1 ) = ^dxexp^-y-r 1 x^ (12) 

that correspond to two basically different choices of the integration contour C in Q. The contour has 
to be chosen so that the integrand vanishes at its ends, i.e. the contour should go to infinity where 
Sftex 3 < 0. One of the possibilities is to choose the imaginary axis as C. This gives the standard Airy 
function 



DC 



X 3 



Ai(Ti) = / cos I — - Tlx I dx (13) 
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Another independent solution to the Airy equation should be associated with a contour that connects 
other infinities. Say, one can choose the contour that goes along the imaginary axis from +00 to zero 
and then goes along the positive ray of the real axis. We call the corresponding function Ai2(Ti). 2 
Then, the "matrix" model (general) solution is, in our case, 

Z(T 1 ) = eAi(Ti) + CAi 2 (T 1 ) (14) 

where £ and £ are arbitrary constants. Now one can use equations (J5J to generate, recursively and 
unambiguously, Z as a power series in t^s. 

Thus, we can define the "matrix" model partition function Z({tk}) as a solution to the defining 
set of equations, ©. Then, the entire freedom that we have in our "matrix" model is due to the 
constants £, f. Choosing these constants, that is, choosing a formal sum of the distinct integration 
contours C, fixes the "matrix" model solution uniquely. 

2.3 "Matrix" integral 

In contrast to the "matrix" model, the "matrix" integral is defined by an integral. The freedom one 
then has is in choosing the integration contour. Therefore, one can take a basis in the space of all 
solutions by choosing some basis contours. This gives us a basis of "matrix" integrals, or the "matrix" 
integral provided with an index associated with the set of basis contours. Then, constructing linear 
combinations of these integrals, one obtains the general solution to the "matrix" model. 

To perform effective calculations, one has to make a clever choice of the basis of integration 
contours. They should be naturally associated with the asymptotical expansion of the integrals. 
Say, dealing with cubic potential, one would better choose the basis contours corresponding to the 
solutions of the Airy equation controlled by different quasiclassical expansions. Indeed, let us make 
an asymptotic expansion of at large T\. Then, the saddle point equation has the two solutions 
x = ±vTT- Depending on the choice of the integration contour C, one should choose one or the other 
solution and expand the integral around this solution to obtain the asymptotic expansion. 

Generally, there are p — 1 solution to the saddle point equation V'(x) = 0, and exactly that many 
different quasiclassical expansions and basis contours. 

Now, in order to obtain Z as a function of all tf~, one may use two different strategies. First of all, 
one can just iteratively solve equations (J5J). The other possibility is to calculate moments @ using 
properties of the corresponding "matrix" integrals, i.e. in the cubic case, those of the Airy functions. 

There is also a more tricky possibility which is in no way technically simpler, but will be of great 
use later. That is, let us consider a possibility of immediate calculating integral This integral is 
typically not a power series in the higher couplings. Indeed, since the degree of potential is now much 
higher than p, there are much more possibilities of choosing the integration contours (=quasiclassical 
regimes, =solutions to the saddle point equation). We have to fix the integration contours so that to 
have a smooth (power) behaviour upon bringing higher couplings to zero, while the behaviour of the 
integral w.r.t. the first p couplings remains arbitrary. This leaves us with the freedom of exactly p 
different integration contours. 

Let us see how it works in the cubic case. Here we have two possibilities of choosing the integration 
contours. Note that the Airy equation (jTT|) can be reduced to the Bessel equation and, therefore, 
its solution is expressed via the cylindric functions. More concretely, the function y /7 T7-Zi/3(§^~f^ 2 ) 
solves the Airy equation, |15j . where Z 1 / 3 (z) is any cylindric function of order 1/3. Now let the 
two basis cylindric functions be the Hankel functions of the first and second kinds H^' 2 \z) (this 
choice corresponds not to Ai(Ti) or Ai2(7i), but to their linear combinations). We also restore the 

2 The difference Ai 2 (Ti) - Ai(Ti) is usually denoted Bi(Ti), see [H>], formula (10.4.33). 
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dependence on all the three couplings (and slightly rescale them for the sake of convenience). Then, 

Z(T 1 ,T 2 ,T 3 ) = J dxe W {iT 3 x 3 -T 2 x 2 + iT lX ) = ^jZ 1/3 (^ 2 ), r, = — - — |^ (15) 

and, using the asymptotic expansion of the Hankel functions, we find that the smooth behaviour under 
T3 — * 0, i.e. 77 — > 00 is celebrated with H@\z) (see ^H], formula (7.13.2)), 

r(2 ),^ ./~2~_/' Az-2tw-it\ ^ (§,m) 1 _ r (§ + "») 



^ ( Z ) ~ J — exp(-i— — — - Y j^r-, (-,m) = — ^ ^ (16) 

i/3^™V^ p v 4 ;^ (2«o-' l 3' j m ir(|-i 1 



Therefore, we finally have (here we put T\ = for the sake of simplicity) 

1 V ■v / \ tvi { 1 \ / 2 7T^: 

m=0 



2 X m 



i.e., with this choice of solution, Z is indeed a power series in T3. One can easily check that the 
corresponding power coefficients coincide with the corresponding moments of the Gaussian integral. 
Therefore, one may really calculate moments of the Gaussian integral in this tricky way. 



2.4 Resolvent and loop equation 

To conclude our discussion of the toy example, we comment on the properties of the resolvent (jlUj) . 
We start from the simplest Gaussian potential. Then, the resolvent 

00 e — T2X 2 j-oo e — T2X 2 poo e~ x2 



G(z)= dx = 2z dx— - = w {y / T 2 z), w(z) = dx (18) 

J — 00 Z X J— 00 Z X J —00 z X 

can be calculated in two different ways. First of all, one can use the formula 

1 



due -u{z-x) ( 19 ) 



Z- X Jo 

and further calculating the Gaussian integral to express the resolvent through the error function 3 , 



w(z) = e z2 



f z 2 \ 2 f°° _ 2 

-iw + 2^/ti J e x dxj = e 2 erfc(-iz), erfc(z) = -i2^ J e x dx (20) 



The other possibility is to find the differential equation for w(z). To this end, let us again use vanishing 
of the integral of the full derivative, 



d ( e' x2 



= / dx 7T ( 21 ) 

J -00 ox \z — X J 

This leads to the equation 

w'(z) = 20F- 2zw(z) (22) 

Solving this equation, one can arrive at l)20|) again. 

In the case of more complicated potentials V(x), one can not manifestly calculate the integral for 
G(z). Moreover, the differential equation is also quite complicated, with order increasing as the degree 
of the polynomial V{x) increases. However, there is a universal form of the equation for the resolvent 
that we discuss now. To obtain it, let us consider the next simple example of the cubic potential. 
Then, we can write, as above, 

(4=^1 <23> 



3 One should be careful with integrating around the point z — x. 
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which is equivalent to 

G'(z) = V'(z)G(z) - [ dxxe^- TlX -4z j dxe^~ Tlx (24) 
Jc Jc 

3 

In contrast to the Gaussian integral, one can not manifestly calculate moments of e~~ TlX . Moreover, 
these moments contain an ambiguity related to the choice of the integration contour. 

Note, however, that the uncalculable part is a (linear) polynomial. Therefore, one obtains that 
the following equation is correct 

G\z) = [V (*)<?(*)]_ (25) 

where [...]- denotes the projector onto negative powers of z, and we took into account that [G'(z)]_ = 
G'(z). Looking at 1)22(1 . one observes that it satisfies (|25j) too. Moreover, repeating the derivation for 
the general polynomial potential, one comes to the same universal result (J23|) . 



3 N x N matrix case 

Now we consider the true matrix integral following mainly the line discussed in the previous section. 
3.1 Loop equations 

We start now with integral (^Q), again not specifying integration contours. (Since this is the multiple 
integral, we have freedom to choose different contour for each integration.) Let us first obtain the set 
of equations satisfied by this integral. As before, we consider the integral of the total derivative 



DMtr 



M n+l exp tr V(M) + V t k ti M k 



dM t i - - ' ^ 

V \ k 



(26) 



where M is the transpond matrix. These equations lead to the following set of constraints (=Schwinger- 
Dyson equations, = Ward identities) |12j 

d d 2 

L m Z N (t) = 0, m>-lL m = Y^k{t k )- + l^r ( 27 ) 

fc > dt k+m a+b ^ m dt a dt b 

a,b>0 

where ^ = NZ N . Note that this time we can not replace derivatives w.r.t. higher couplings with 
higher order derivatives w.r.t. the first coupling. As before, in order to have non-trivial solutions to 
these equations, we need to shift first p couplings: t k — > T k + t k , k = l,...,p and then deal with the 
partition function Q} asa power series in the couplings t k s. Then, one also needs to add to constraints 
(|27j) the conditions 

W L = § L Vfc = 0,...,n + 1 (28) 
oT k dt k 

However, since now we have the set of equations w.r.t. various couplings, they can no longer be 
reduced to a single equation w.r.t., say, T\. Therefore, the freedom in solving equations (|2*T|) is much 
larger. 



3.2 Matrix model 

How many solutions to equations (|27j) do we expect now? To understand this, let us again study 
an oversimplified example. Namely, look at the model similar to P|. but without the Vandermonde 
determinant. Then, it reduces to the product of independent factors 
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Each factor here is a solution to the corresponding ordinary differential equation. Say, in the cubic 
case, one has 

N 

Z N (T\t) = J] [fcAifa) + QM 2 ( V )] = ]T ~ fe [Ai( V )] k [M 2 (f})] N ~ k (30) 

i k 

where 3^ are the coefficients constructed from products of and Q's which determine the freedom 
in the matrix model partition function. Since are arbitrary coefficients, one may interpret them as 
counterparts of Fourier series coefficients of an arbitrary function Zn(ij) of ij, the Fourier exponentials 
being substituted with the combination of the Airy functions. 

Generally, as we discussed in the previous section, for the polynomial potential of degree p, there 
are p — 1 independent basis functions, i.e. the freedom in the matrix model partition function is an 
arbitrary function of p — 1 variables. The same counting certainly remains valid for the model with 
the Vandermonde determinant present, 

Let us now understand the origin of the arbitrary function of p — 1 variables in terms of constraints 
(|27[). To this end, note that if one solves them recurrently, expanding the partition function into the 
power series in couplings t^s (like (|7jl). only the first two constraints of 1)27)1 are really restrictive for 
Zn(T), while all other constraints just allow one to restore recurrently the dependence on the couplings 
i/c's. As we already explained, (|27|) are less restrictive than © since they can not be reduced to an 
ordinary differential equation in T\. 

The first two constraints of 1)27)1 are linear in derivatives and, therefore, we consistently truncate 
them to t = and then express two derivatives, say, dZ/dT p+ \ and dZ/dT n through dZ/dTi with 
I = 0, . . . , p — 1. As a corollary, the partition function can be represented as 

Z N (T) = f dkz(k; rj p )e^ kx - k2 ^ (31) 



with an arbitrary function z of p arguments (k, rj2, ■ ■ ■ , r] p ). Here the following variables invariant w.r.t. 
the first constraint are used, 

1 . . 

w = ^-py log T p+ i, x = T + . . . ~ rjp + i (32) 

r) , _ ( T k | k ( k ~ 2 ) ! v-y ) i (p+ 1 ) l (p- l V- T T k-i-\ T i \ T -jfi ( oo\ 

m - \ +p + — -\ — (k-l-l)\ p p p+1 J p+l 1 ' 

We discussed a particular case of these formulas in s.2.3. 



3.3 Matrix integral 

Let us discuss now what could be a choice of basis matrix integrals. We basically repeat the procedure 
we applied in the ordinary integral case, i.e. use the saddle point approximation. Different saddle 
points M = Mq are given by the equation V'(Mq) = 0. If the polynomial 

p 

V'(x) = l[(x-a i ) (34) 
i=l 

has roots aj, then, since Mq are matrices defined modulo [/(iV)-conjugations (which allow one to 
diagonalize any matrix and permute its eigenvalues), the different saddle points are represented by 

Mo = diag(ai, . . . ,ai;a>2, ■ ■ ■ ,ot2', ■■■ ;a p ,...,a p ) (35) 

with a.i appearing iVj times, X^=i ^ = ^ • f &c t> there is no need to keep these iVj non-negative 
integers: in final expressions they can be replaced by any complex numbers. Moreover, iVj can depend 
on Tfc (i.e. on the shape of V(M)). 
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Now, using at the intermediate stage the eigenvalue representation of matrix integrals, one can 
rewrite 0J [7] the matrix integral (??) over N x N matrix M as a p-matrix integral over Ni x iV, 
matrices Mj (each obtained with the shift by a^: just changing variables in the matrix integral (??)), 
which is nothing but the multi-cut solutions [2] 



Z v (t\M ) 



p 

n 

i=l 



n™,exp(^tr4 l) MMn« 



2NiNj 
ij 



i,k 



i<j 



«P ( 2 f (-) fc (fe ti,,? ! t^Mftr J Mj 



(i) 

The variables t k are given by the relation 



V / oo 



fc=0 \i=l / i=l \fe=0 



(36) 



(37) 



with arbitrary JVj x iVj matrices Mj. 

Thus, we finally have p — 1 basis functions that can be described by the proper choices of the 
integration contours, and are associated with different solutions to the saddle point equation. 



3.4 Resolvent and loop equation 

Another form of constraints (|27[) is produced by rewriting the infinite set through the unique j 
function of all single trace correlators 



k>Q 



_l d_ 



T = g 2 \ogZ v 



generating 
(38) 



Introducing the notation v(z) for J2k^k zk : one obtains the loop equation 

'v'(z)pM(z\t)] = V'{z)p i - l \z\t) - \v'(z)p w (z\t) 



(39) 



(p«(*|t) 



v'{z)p {1 \z\t) + g 2 V(z)p^(z\t) 



In order to consider (connected) multi-trace correlators, one needs to introduce higher generating 
functions (also named loop mean, resolvent etc) 



Am) 



(zi, z m \t) = V(zi)...V(z m )J r 



Note that G(z) introduced in (|1U|) is equal to Z(T\t)p^(z\t). Moreover, taken at all tj. 
reduces to (|25|) . However, the quantity p^\z\0) generates only correlators of single-trace 
(moments), which is not enough in the matrix case. Therefore, in this case one needs to 
whole quantity t). It can be, however, expressed through p^ m \zi, z m \0), 

pW(z\t) = "[/•••/ V ^ ■ ■ ■ v(Zm)p {m+1) (z, Z X ,..., Z m \0) 

m>0 



(40) 

= o, (ESJ 

operators 
know the 



(41) 



4 Genus zero solution 

One can solve the loop equations (|3*9*|) recurrently (see, e.g., [S]) expanding them intosum of p( m \zi, ...,z m \0) 
and these latter into series in g. This gives one the double recurrent relation (in m and order of g). 
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Note that the solution of these recurrent relations in the leading order can be immediately obtained. 
Indeed, omitting the last term from (|39j) . one obtains that 

p m) {z) = v '(z)-y(z) (42) 

(where the first superscript refers to the leading in g approximation) with 

y\z) = {V'{z)f-AP p ^(z) (43) 

where P p -\(z) = V (z)p^ (z\0) is a polynomial of degree p—1. Coefficients of this polynomial 
depend on T^'s and first derivatives of the arbitrary function of p — 1 variables, which parameterize 
solutions to constraints (|2T)) (=the loop equations). 

Note that formula (|l3*|) gives a hyperelliptic Riemann surface of genus p—1, since mapping gives 
p cuts on the complex plane (this is the celebrated multi-cut solution we discussed in the introduction) . 
It turns out that, at least with some specific fixing of the ambiguity in solutions of the loop equations, 
the (multi-)resolvents can be constructed as differentials given on this Riemann surface, |17U14j . This 
specific fixing is given by the conditions 

p( m+1 \z, Zl ,...,z m \0)dz = (44) 

Ai 

for all the (multi-)resolvents except p^\z), 

p^ 1 \z) = S i (45) 

Ai 

where Si are arbitrary constants not depending on and Ai are ^-cycles on the Riemann surface 
(|43j) . The constants Si can be, in fact, associated (up to the factor g) with N% from (|36j) . Moreover, 
the basis functions (|3*B]) exactly lead to conditions (|44 |) -(|45 |l . These conditions are also associated with 
Seiberg-Witten-Whitham system corresponding to the Riemann surface (|43|) |13j . 

We do not go into further details here, just referring to ^U] and ^1] for the latest development 
and proper references. Let us only note that conditions (|44 |l -(|45 |l distinguish specific solutions that 
survive while smooth changing the number p of couplings T/., much similar to what we observed in 
s.2.3. Indeed, for any number q of non-zero TVs, one may require the curve to be of the form 

y\z) = {V'{z)f - 4P,_i(*) = H 2 q _ p (z)R 2p (z) (46) 

where H 2 _ p (z) and R2 P (z) are polynomials. Therefore, independently of q, one has the hyperelliptic 
curve y 2 (z) = R2p(z) of genus p—1. Moreover, the freedom one has in matrix model in this case is 
dictated by p — 1 constants Si, see the details in |14| . The crucial difference from the case considered 
in s.2.3 is that, when one is dealing with the ordinary integral, one may calculate higher moments 
differentiating the integral in lower couplings several times, while, in the matrix case, it is not possible 
because of the traces. Therefore, technically this is very convenient to add traces of matrices in higher 
degrees directly into the matrix potential. One, however, has to be careful not to change the solution 
of the matrix model with this procedure. This is what exactly achieved with preserving the curve 
y(z), ©. 
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